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We investigate the correlators in unitary minimal conformal models cou- 
pled to two-dimensional gravity from the two-matrix model. We show that 
simple fusion rules for all of the scaling operators exist. We demonstrate the 
role played by the boundary operators and discuss its connection to how loops 
touch each other. 
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1 Introduction 



The understanding of two-dimensional quantum gravity has experienced great 
progress through the study of the matrix modelsB The one-matrix model has infinite 
number of critical points which are considered to represent the {2m + 1, 2) minimal 
conformal models coupled to two-dimensional gravity. The two-matrix model has 
critical points which correspond to the (m + 1, m) unitary minimal conformal models 
P, 0, ^. In this paper we investigate the unitary minimal model (m + 1, m) coupled 
to two-dimensional gravity from the two-matrix model. 

The emergence of the infinite number of scaling operators aj is one of the most 
important properties of the matrix models. Before coupled to gravity, the minimal 
model has finite number of primary fields. After gravitational dressing, however, 
infinite number of scaling operators emerge. This phenomenon can be understood 
as follows. In the Kac table we can divide the primary fields $r,s into those which 
are inside the the minimal conformal grid 1 < r < q — 1, 1 < s < p — 1 and those 
outside, which correspond to the null states. Before dressed by gravity, the fields 
outside the minimal conformal grid decouple from physical correlators. After 
gravitational dressing, they cease to decouple [|], and become infinite number of 
scaling operators. The similar phenomenon has been shown in continuum frame- 
work. Through the examination of the BRST cohomology of the minimal model 
coupled to Liouville theory, infinite physical states were shown to exist [§. These 
states have their counterparts in the matrix models as the scaling operators. Some 
of the scaling operators do not have their counterparts in the BRST cohomology, 
which we will discuss later. 

In ordinary {p, q) minimal conformal model the primary fields satisfy certain 
fusion rules p; three-point function {^ri,si^r2,s2^r3,s3) is no n- vanishing only when 

1 + |ri — r2| < T3 < mimjri + r2 — l,p}, ri + r2 + r^ = odd 
1 -h |si - S2I < S3 < mim{si + S2 - I, g}, Si -F S2 + S3 = odd . (1.1) 

It is interesting to examine how the fusion rules change when the matter couples 
to gravity. The three-point functions involving lower dimensional scaling operators 
were examined from the point of view of the generalized KdV flow in in the 
case of {m + 1, m) unitary matter. It was shown that the gravitational primaries 
aj (j = 1, ■ ■ ■ , m — 1) satisfy fusion rules of BPZ type; for ji + J2 + js < 2in — 1, 
(cTji crj2 C"j3 ) is non- vanishing only when 

I + IJ1-J2I <J3<J1+J2-1. (1.2) 

The fusion rules were also examined in continuum framework |]^. As for the grav- 
itational descendants, however, we think clear results have not been obtained. In 



^See for example H] for review. 



this paper we would like to clarify the fusion rules for all of the scaling operators 
including the gravitational descendants in the case of unitary minimal model. 

Macroscopic loop correlators, which are the amplitudes of the surfaces with 
boundaries (loops) of fixed lengths, are the fundamental amplitudes of the ma- 



trix models. It was shown [10, 11 1 that these correlators have more information 
than those of local operators and that the latter correlators can be extracted from 
the former correlators explicitly in the case of c = 0, 1/2, 1. They argued there 
that macroscopic loops could be replaced by a sum of local operators whose wave 
functions satisfy the Wheeler-DeWitt equations. 



We generalize this idea to the loop correlators g, |T^, |1^, Q in the cases of the 
general unitary minimal models, and derive the fusion rules for all of the scaling 
operators. First we derive the explicit form of the expansion of loops in terms of 
the scaling operators, and then deduce the three-point correlators from the loop 
correlators which were calculated in [|T3l HM from the two-matrix model. We show 



that the three-point correlators of all of the scaling operators satisfy certain simple 
fusion rules and these fusion rules are summarized in a compact form as the fusion 



rules for three-loop correlators |13 



In matrix models, there are infinite subset of the scaling operators Bj {j = mod 
m+1) which do not have their counterparts in the BRST cohomology of the minimal 
model coupled to Liouville theory. In the case of one-matrix model, Martinec, Moore 
and Seiberg [^ argued that these operators are boundary operators which couple 
to the boundaries of two-dimensional surface. They proved that one of them is 
in fact a boundary operator which measures the total loop length. But little has 
been discussed on the role played by the rest of these operators. We examine the 
geometrical meaning of these operators and its connection to the touching of loops 
in the case of general unitary models. We show that the boundary operators have 
the role to connect several parts of the loops together. We also discuss the relation 
of the boundary operators to the Schwinger-Dyson equations proposed in |T^ . 



2 Expansion of loops in local operators 



We consider the {m + l,m) unitary minimal model coupled to two-dimensional 
gravity from the two-matrix model with symmetric potential. The partition function 
Z is defined by 

e^ = |did5e-xtK^(^)+f^(^M^) , (2.1) 

where A and B are hermitian matrices, and f/ is a certain polynomial. In this article, 
we limit our discussion to the two-matrix model with symmetric potential and to 
the critical points which correspond to the unitary minimal models. In the case 



of asymmetric potential, some of the boundaries (loops) of two-dimensional surface 
would have fractal dimensions different from the usual dimension of length. 

In [0, in the case of the one-matrix model, it was shown that the loop operator 



can be expanded in terms of local operators, that is, the loop can be replaced 
with the infinite combination of local operators, except some special cases. It was 
argued that this is the case for the general {m + 1, m) unitary models. Whether this 
replacement can be done safely or not is connected with whether the corresponding 
classical solution exits or not in the limit of small length of the corresponding loop. 
This claim is quite natural because all of the scaling operators are expressed in term 
of one matrix A as aj = Tr(l — A)^^^!'^ = X^n'^nO)^^^ Tr A"' in the one-matrix 
model. 

In the two-matrix model, it appears that this idea is not the case since the 
direct connection of the scaling operators to the operators Tr A" or Tr i?" is not 
clear. But this expansion is considered to be possible by the following reason. When 
one of the loops on two-dimensional surface shrunk to a microscopic loop, the loop 
represents local deformation of the surface. The microscopic loop can be replaced 
by the insertions of local operators. The loop correlators except one-loop case are 
continuous when the length of one of the loops approaches zero, so that we expect 
that a macroscopic loop can also be replaced by a sum of local operators. 

In this section we derive explicitly the expansion of loops in local operators in 
the case of the unitary minimal models. Using this relation, we can deduce the 
amplitudes of local operators and those involving both loops and local operators 
from the loop correlators. 

Let us recall that the two- loop correlators in the (m + l,m) unitary minimal 
model coupled to two-dimensional gravity are [^ |1^ 



1 M 



m 2 ii + (±)"^£2 



1^2 Y. {±f-^K^{Mh) K,_^{Mi2) . (2.2) 



k=l 



for £i ^ £2- Here w~^{£) and w {£) represent loop operators which create loops made 
by the matrices A and B respectively, and we introduced a notation 



sin7r|p| 
7r/2 



K,iM£) = j^K,{M£) . (2.3) 



Here the parameter M is defined by 



-tt) =— XT' A-A* = -«V 2.4 

2 J m + 1 

where A* represents the critical value of the bare cosmological constant, and fi is 
the renormalized cosmological constant. 



The definition of a set of the scaling operators has arbitrariness which conies from 



the contact terms [|10|- As a set of local operators, we take the scaling operators dj 
whose wave functions \I/j(^) satisfy the (minisuperspace) Wheeler-DeWitt equations 



'§^P'^''<£^ 



M/,(£) = . 



(2.5) 



It was shown flQ] that these scaling operators correspond to the dressed primary 
fields of the conformal field theory in the case of one-matrix model. In terms of these 
scaling operators, we can obtain simple fusion rules for three-point correlators, which 
we will show in the next section. 

We normalize the wave function of dj as 



(a, w^{i) 



7 /M\™ ~ 

— — KAMtj , J>1 ^O(modm) . (2.6) 

m \ 2 J m 



Note that the normalization factor sin -^vr in Kj{Mi) in eq. (|2.6|) is consistent be- 
cause there are no scaling operators dj for j = (mod m) in the matrix model. We 
can express the right hand side of eq. (|2.2| ) as an infinite sum in terms of Kj{Mi,2) 
for li < (,2 (see appendix A): 



w;+(£i)w;±(4) 



m—l oo 



m 



E E (± 

fc=l n=— CO 



>fc-l 



+ 2n 



m 



Lk 



+2n 



(M£i) Kj,.,jMi2 



(2.7) 



Comparing eq. ( p.6|) with eq. ( p.7|) , we expect the following expansions of the loop 
operators in term of the local operators: 



" -^W = - E E i±r' (f) 'i^ 

m f-^ „_ _ V 2 / '™ 



+2n\ 



(Mi) ai 



k+2mn\ 



(2.8) 



fc=l n=~QO 

These expansions are the generalizations of those in the case of the one-matrix model 



10| to the cases of general unitary minimal models coupled to two-dimensional 



gravity. 

Since the loop correlators are symmetric under the interchange of two kinds of 
loops, that is, {w^{(ii) ^^(£2)) = (""^"(^i) ""^"(^2)), the wave functions of the scaling 
operators with respect to loop w^(£) are read as 



(a, w-{l) 



^y-' {a, w^ii) 



(2.9) 



3 Fusion rules for scaling operators 



Using the expansion of loops eq. (|2.8| ), we can obtain the correlators of the 
scaling operators from loop correlators.^ In this section, we show that there are 
rather simple fusion rules for all of the scaling operators. The fusion rules involving 
the gravitational descendants {ctj, j > m + 2 ) have not been clear from the point 
of view of generalized KdV flow. 

3.1 One- and two-point functions 

Let us examine one- and two-point functions first. Since the one- loop ampli- 
tude diverges when the loop length approaches to zero, this amplitude include the 
contribution which is not represented by the local operators. Extracting the parts 
proportional to Ii, {u > 0), which parts can be considered as the contributions from 
local operators, from the one-loop amplitude 



w^{i) 



} = 


f 

1 + 
V 


mj 


t 


-1 (M-^ 
\ 2 y 


1^1 


2 , 


2+i 


[h. 




(M£)- 


-1-2 



.{Mi) 

_{Mi) -U{Ml) 



.{Mi) 



we can obtain the one-point functions of the scaling operators 

{aj) = , J ^1,1 + 2m. 



(?i) 



Let us turn to the two-point functions. Substituting eq. 
obtain the two-point functions 

.^j. 2j/m 

(a^dj) = 6ij j ( — j , «, j 7^ (mod m) 

Note that we obtain diagonal two-point functions. 



into eq. 



(3.1) 

(3.2) 

(3.3) 

15[), we 

(3.4) 



3.2 Three-point functions 



The three-loop correlator from two-matrix model is obtained in |jT3 

'M^ 1-^ 



as 



W+{i^)w+{i2)w+{is) 



1 



■m{m + 1) 



'It2t3 



^The multi-loop correlators were examined in ||lj] from two-matrix model. These correlators 
were also examined in [ha from the viewpoint of random surfaces immersed in Dynkin diagrams. 



X 



J2 K^_^{Mi,) K^_^{Mi2) K^_._,iMi3) , (3.5) 



(fcl-l,fc2-l.'=3-l) 



Here we have denoted by V^' 



Oi — aj > for j = 1 ~ 3 , 



^3 = \ (01,02,03) I J2 

^ aj = even < 2(m — 2) , Oj = 0, 1, 2, ■ ■ ■ > 
i=i J 



(3.6) 



It was shown [|l^] that the selection rules in eqs. ( p.5|) and (|3.6| ) correspond to the 
fusion rules [|, for the dressed primaries ( 0ij, i = 1, ■ ■ ■ , m — 1) by studying 
the small length behavior of the three-loop correlator. Using the expansion of loop 
eq. ( ^.81 ), we can show that the selection rules in eq. ( |3.5| ) represent the fusion rules 



for all of the scaling operators in a compact form. Let us show this in the following. 
Using the formula, 

°° \p + 2n\ 

zKi_\p\{z)=7l Y. 1 -l^J\p+2n\iz) , (3.7) 

^r^^ smnlp + 2n\ 
we first expand the three-loop correlator eq. (|3.5|) as 



'm[m-\-i) \ Z / „(m) ni=-oon2=-oon3-oo 

(^ + 2^1) (t + 2r^2) (^ + 2n3)/|^^,„^|(M£i)/|^^,„^|(M^2)/|M+2n3|(^^3) • 

(3.8) 

Comparing eq. (|3^ ) with eq. (|2.8| ), we can extract the three-point functions 

{'^\ki+2mni\'^\k2+2mn2\'^\k3+2mn3\) 

= Ck,k2ks^-—TTl[iki + ^mrii) — , (3.9) 

m{m + l) ^Ji \ 2 / 



where 



^'^'^'■' -jo, otherwise " ^'^'^^^ 



For rii = 0, eq. ( p.9| ) is nothing but the correlator of the gravitational primaries. 
For the gravitational primaries, eq. (|3.4|) and eq. ( p.9|) agree with the correlators 



obtained in |^ from the generalized KdV flow up to a factor —2. Note that we obtain, 
here, the correlators of the gravitational descendants as well. In 0, the fusion 
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rules for the gravitational primaries were examined in continuum framework. We 
have found here the fusion rules for the gravitational descendants as well as for the 
gravitational primaries. These fusion rules are similar to those for the gravitational 
primaries due to the factor Ckik2k3 in eq. (|3.9|). 

Introducing the equivalence classes [dk] by the equivalence relation 

dk ~ ^\k+2mn\ , n e 2 , (3.11) 



we can consider the fusion rules in eq. (|3.9| ) as fusion rules among [dk] {k = 1, ■ ■ ■ , m— 
1). Note here that the class \ak\ does not correspond to the set which consist of the 
gravitational primary Ok and its gravitational descendants aiiOk) , / = 1, 2, ■ ■ • in 
P introduced from the viewpoint of KdV flow. 



3.3 Further on the fusion rules 

In this subsection, let us examine the fusion rules in eq. ( p.9| ) further and summa- 
rize the relation of the scaling operators to the primary fields in the corresponding 
conformal field theory. 

In the (p, q) minimal conformal model, the primary field $rs has the conformal 
dimension 

^^^=^^-^^f-^P-^? , (3.12) 

Apq 

where r and s are positive integers. Since we have A,. ,, = A^+g ,,+p = ^q-r,p-s^ the 
corresponding primary fields can be regarded as the same one. The integers r and 
s can thus be restricted in the range 

l<r<g — 1, l<s 7^0 (modp), pr < qs (3.13) 

(see fig. |I] ). In fig. [^, the primary fields in the region ((2)) or ((2))' are the secondary 
fields of those in the region ((1)). In general, the fields in the region ((n+l)) or ((n+1))' 
are the secondaries of the fields in ((n)) or ((n))'. Since the secondary fields correspond 
to null vectors, those fields decouple. One can thus construct consistent conformal 
field theory which include only the primary fields in the region ((1)) (i.e. inside the 
minimal table), that is, the (p, g) minimal model [^]. Coupled to Liouville theory, 
however, the fields outside the the minimal table fail to decouple and infinite 
physical states emerge accordingly. These states are considered to correspond to 
the primaries outside the minimal table. This correspondence is implied by the 
BRST cohomology [§ of the coupled system. 

In the minimal model coupled to Liouville theory, the gravitational dimension 
A^^ of the dressed operator for (^r,s is given by the relations 

aG ^ 1 _ ^ »r.. _ p + q-\pr-qs\ 
7 7 2g 



q 




2p 



3p 



Figure 1: the range of (r, s) 



where r and s take the values in the range eq. (|3.13| ). On the other hand, in the 
matrix model, the corresponding relation for the scaling operator dj is given by 



Oj ^ p + q-j 

7 2g 



From eq. ( |3.14| ) and eq. ( ^.15| ), we should take as 

j = \pr -qs\ , j = 1, 2, ■ ■ ■ 7^ (mod q) 



(3.15) 



(3.16) 



for Bj. 

Consider now the relation of a\k+2mn\ to the primary field $r,s of the unitary 
{m + l,m) minimal model. Let us first compare the two sets 



Sk = \\k + 2nm \ | n G 2 > 



(3.17) 



and 



I |pr — gs| [ = I |(m + l)r — ms| > 
= I r' + (s — r — l)m > 

where r and s are positive integers in the range 

l<r<m — 1, l<s, r + 1 < s 



(3.18) 



(3.19) 



and r' = m — r. Note that we include s = (mod m + 1) here. Decomposing the 
set Sk into two sets as 

Sk = S+®S^ , (3.20) 



O; O2 Oi 



o 
® 



0"|l+2mn| 
<7|2+2mn| 
C"|3+2mn| 



® 




10 



15 



HO,) 



Figure 2: the scaling operators d\k+2fnn\ in the (5, 4) minimal model coupled to 2d 
gravity 



where 



St 



\k + 2nm \ n = 0,1,2, ■■■ } , 

{(m- A;) + (2n' + l)m | n' = 0, 1, 2, ■ ■ ■ } , 



(3.21) 



and comparing eq. ( p.21|) and eq. ( p.l8|) , we can express the sets S^ and S^^ in terms 
of |(?7i + l)r — ms\ as 



Sk 

Oh. 



I \{m + l)r — ms\ \ r' = A;, s — r = 2r2 + 1, ra = 0, 1, 
I |(?Ti + l)r — ms] I r = k, s — r = 2n' + 2, n' = 0, 1, 



}. (3.22) 



From eq. (|3.22| ), the following correspondence is obtained: 



H+2mn\ {n > 0) 
^\k+2m{-l~n')\ {n > 0) 



^ ^m^k, r+2n+l ijl > 0) 
^ $fc, r+2n'+2 (^' > 0) , 



(3.23) 



where s 7^ (mod m + 1). As for the scaling operators corresponding to s = (mod 
7T7, + 1), we will discuss these in the next section. 

As an example, we depicted the scaling operators on the r-s plane for the case 
of m = 4 in fig. ^ In this figure we showed the equivalence classes [a^] explicitly. 



4 Boundary operators 

4.1 Boundary operators and touching of loops 

The scaling operators dj {j = mod m + 1), which correspond to s = (mod 
??7, + 1) on the r-s plane do not have their counterparts in the BRST cohomology 



of the minimal model coupled to Liouville theory. In |T^ it was proposed that the 
scaling operators which do not occur in the BRST cohomology of Liouville theory 
are boundary operators and one of them, which is a^ = ai(Ci) in the case of pure 
gravity, was in fact proven to be a boundary operator for the one-matrix model and 
the Ising model case. We would like to examine the role played by the operators 
d'n{m+i), n = 1,2,- ■■ 7^ (mod m) as well as Bm+i for general unitary minimal 
models. 

Let us denote these operators by 

Bn = an{m+i), u = 1,2, ■ ■ ■ ^ {mod m) . (4.1) 

In the matrix models the loop amplitudes contain the contribution from the surfaces 
with loops touching each other. In two-loop case, let us consider the surfaces in 
which the two loops touch each other on n points. When we shrink one of the 
loops to a microscopic loop, the other loop splits into n loops, which are stuck each 
other through the microscopic loop (see figs.H, ^ and ^. Since the microscopic 
loop represents a sum of the scaling operators, the wave functions of some scaling 
operators contain the contribution from the surfaces with split loop. 

We now show that the boundary operators indeed represent these surfaces. From 
eqs. ( [^.6D and (|3.1|), the wave function of S„ and the one-loop amplitude are 

B^w^i)) = nil + -) (-) "^ K„(,+i)(M£), (4.2) 

(w^i)) = (1 + ^) r^ (f ) '^"^ K.^^m) ■ (4.3) 

We denote the Laplace transformation of any function f{ii,i2, ■ ■ ■) oi loop lengths 

by 

C [fih, i2, ■■■)]=/ di, / di2--- e-(^^f^+^^^^+-)/(£i, £2, ■■■). (4.4) 

Jo Jo 

In the space of Laplace transformed coordinate (, we have 

i-\BnW^{i))\=-(—j 2coshn(m + l)^, (4.5) 

C[{w+{i))\=-i — ] 2cosh(m+l)^, (4.6) 

10 



where we have used the relation 



£ 



-rV|ir^(M£) 



2 cosh mi>6 



(4.7) 



and ( is parametrized as ( = McoshmO. Note here that w^{i) represents a loop 
with a marked point and i~^w~^{i) represents a loop without a marked point. Since 
cosh n(m + 1)6' can be expressed as a polynomial of cosh(m + 1)^, 



2coshn(?7i + 1)6' = 2 T„(cosh(?7i + 1)6' 



[n/2] 

Y, 4"^ [2 cosh(m + 1)6 

r=0 



"-2'- („) _ (-l)""^ fn-r 
n — r \ r 



{A.i 



where T„ is the Chebeyshev polynomial, we obtain the following relation: 



C 



r=0 



2 / 



{c[-{w^i))]}' 



n-2r 



(4.9) 



In the space of loop length, the above relation means that the wave function of Bn 
is equivalent to a sum of the convolutions of disk amplitudes: 



Bnw+{e) 



[(n-l)/2] 

' E 

r=0 






V 2 / 



n— 2r 



(_l)"-2r 0^+ (£) _ (410) 



Here we introduced a notation ^j'" = {w^{i)), and 0^i" (£) denotes the convolu- 
tion of s Ai{(iys, for example 



0^+1 {£) 



00 roo 



JO 



d£id^2 S{h +i2- i)At{ii)At{i2) • 



(4.11) 



From eq. (|4.10| ) we can conclude that the operator Bn couple to the point to which 
s {< n) parts of the loop are stuck each other in the case of one- loop amplitudes. 
When there are more than one loop, we infer that the operator couples to the point 
to which s parts of several loops are stuck each other; the operator will not recognize 
that it is touching different loops this time. 



Using the following relation 

[(n-l)/2] 



2 cosh X 
we also obtain 



E 

r=0 



n 



2 cosh(n — 2r)x, (up to constant) 



(4.12) 



QAtY H) = i-ir' E f")(f) {Bn-2rWHi) 



r=0 



(4.13) 
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Here we drop the constant term in eq. (|4. 12|) when we carry out the inverse Laplace 



transformation. From eq. ( [4.13|) , we see that the boundary operator coupled to the 
point on which n parts of loops are touching each other is given by 

[(n-l)/2] ^ 2r(l+i) ^ 

S- = (-!)"« E (")(f) e,._... (4,14) 

Now let us consider the boundary operators when there are two loops on two- 
dimensional surface. As for Bi, we expect that {w'^ {ii)w^ (£2)^1) should be pro- 
portional to {£1 + £2) {w~^{£i)w'^{£2))- Let us confirm this in the following. From 
the three loop correlator ( ^.51 ), the expansion of loop operator ( |2.8D and the wave 
function of a\k+2mn\ ( |2.6| ), we obtain the following correlator with two loops and a 
local operator: 

x£i4(— + 2ns)K^_!niM£^)K^,^{M£2) . (4.15) 

Consider the amplitude for Bi = Bi = Bm+i = 5'|m^_2|. Since Ck-i^^k2,7n-i is nonvan- 
ishing only for the case of ki + k2 = m, we have 

w^{£^)w^{£2)B^ = - Y.{±)'-' { — ) £i£2K^{M£,)K,_^{M£2) . (4.16) 
I m ^ V 2 / 



Comparing eqs. ([4.16|) to ( |2.2| ) we obtain the desired relation 



w+{£r)w^{£2)Br) = K + (-1)™£2} {w+ {£{)w^ {£2)) . (4.17) 

Note here that we have {Biw+{£)) = (-1)™ {Biw-{£)). 

Next, let us consider B2. Since we infer that the insertion of B2 should play the 
role of connecting two parts of loops together, we expect the relation 



w+i£i)w^i£2)B2) = 2£, J^ ' d£[ (w^ i£[)w+ (£2)) {w^i£i - £[] 

+2^2 f' d£'^ {w+{£i)w+{£'^)) {w+{£2 - £'2)) 
+2£^£2(w+{£i+£2)) . (4.18) 



The third term in the right hand side of eq. ( [4. 181 ) represents the contribution from 



the surfaces with loops w^{£i) and ^"'"(£2) touching with each other on a point. 
Let us confirm the relation (|4.18|) in the following. In this case, it is convenient to 
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-^ 




(a) 



(b) 



Figure 3: (a): A surface with two loops touching each other on a point, (b): When 
one of the loops shrinks to a microscopic loop the microscopic loop is equivalent to 
the insertion of the operator denoted by the dot on the loop. 




■^ 




(a) 



(b) 



Figure 4: The case of a surface with two loops touching each other on two points. 




-> 




(a) 



(b) 



Figure 5: The case of a surface with two loops touching each other on three points. 
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consider in the space of Laplace transformed coordinates Q. In this space eq. ( [4.15| ) 
reads as 

{WHC^)WHC2) ?|.3+2^n3l> = ;;^^ (y) (^ + 2^3) 

a_a_f^ . , xfc.^i sinh(m - ki)9i sinh(m - ^2)^2 ) . ^ 

'9Ci<9C2 \fcf^, ^'"'"'^^ sinhm^i sinhm02 j'^ ""^ 

where W'^^Q) = £[w^(£j)]. Due to the relation 



^2 sinh(m — ki)6i sinh(m — A;2)^2 



fci,fc2 



sinh m9i sinh m9 



2 



-1 fsmh{m-k3)9i ^ fe^ sinh(m - fc3)g2 \ ,^r,Q^ 



2(cosh6'i =1= cosh^2) V sinh?7i^i sinh 7/1^2 

we have 

W^{C,)W^{C2) ^lk,^2mn^ = ,, , ,, h^ (- + 2^3) 



2(m + 1) V 2 / 'm 

d d j 1 /sinh(m - /ca)^! y,^ sinh (m - ^3) 6^2 \ 1 

5Ci 5^2 \ cosh 6*1 =F cosh 6*2 \ sinh mOi sinh m6'2 / J 

(4.21) 

Since we should take k^ = 2 for ^2 = —^2{m+i) (for m > 3), we obtain the amplitude 
fori32 



W+iCi)W^iC2)l32] 
Tl fM\^ d d ( 1 /sinh(m-2)ei sinh(m-2)02' 



m \ 2 J dCi dC,2 \ cosh 9i =F cosh 62 \ sinh m^i sinh 17162 / J 

(4.22) 

On the other hand, from the amplitudes H] 



WnC,)W^C2)) = ^4-\n ^°«h^i--°^h^2 



9^1 9(^2 cosh 771^1 — cosh 17102 

d ( 1 sinh6'i 1 



8(2 [ cosh ^1 — cosh 02 mM sinh m0i (^i — (2 
we obtain the relation 

-2 A {(w^+(co#-^(C2)) (w^c,))} + (1 - 2) 
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(4.23) 



-2^ ' 




90 9C2 

2 /M\- 9 9 [ 1 /sinh 6*1 cosh (m + 1^ 



m \ 2 J dCi dC,2 \ cosh 6'i — cosh 62 \ sinh mOi 




One can easily show that the right hand side of eq. ( 4.24 ) agrees with that of 



eq. ( [4.22|) . Putting eqs. (|4.24|) and (|4.22|) together and performing the inverse 



Laplace transformation, we obtain the desired relation eq. ( |4.18| ). 
As for {w'^ {(ii)w^ {(^2)^2) , from the amplitude [Q 

W+{Ci)W-{C2)) = TT-T^ln(cosh^i + coshes) , (4.25) 

Cc,i ac,2 

we obtain the relation 

(w+(£i)w-(£2)S2) = -2£i Ij' M[ {w^{e^)w-{(i2)) {w^{(.i - (■[)) 

+2^2 /'' M'^ (w^{li)w~{Q) (w~{(i2 - i'2)) (4.26) 

in a similar way. In this case, the operator B2 does not connect the different kinds 
of loops w~^{ii) and ^"(^2) together. 

We have shown that the operator B2 connects two parts of the same kind of 
loops together in the case with two loops. We infer that similar phenomena occur 
in general; the operator Bn connects n parts of the same kind of loops together in 
the case with any number of loops. 

4.2 Connection to the Schwinger-Dyson equations 

We can observe a close relationship between the boundary operators and the 
Schwinger-Dyson equations proposed in [|l^. Continuum limit of the Schwinger- 



Dyson equations for loops in the two- and multi-matrix models were proposed in 
17 1 under some assumptions. It was shown [|l^ that these equations for the two- 



matrix model contain W3 constraints, which were derived explicitly in [|I^. The 



integrability of these equations were shown in [|I9|. These facts justify the proposed 
Schwinger-Dyson equations. 

Let us consider the connection of the boundary operators with the Schwinger- 
Dyson equations. For the {m+1, m) minimal models, the following Schwinger-Dyson 
equations were proposed in p 



d£' {w^^\£')w^^^ {i - i'; [n{a)y) w^^^li) ■ ■■w^^\L))' 
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+gY,l, (^« (£ + £,; [n{a)Y) w^^\l{) ■ ■ ■i/;«(^,_i)i/;«(£,+i) ■ ■ ■ u;«(4))' 

i 

for J = 0,---,m-2 , (4.27) 

and 

(m;(^) (£; [Hia)]^-^) w^^\ii) ■ ■ ■ u^^^HC))' ~ . (4.28) 

Here (• ■ ■) represent loop correlators that are not necessarily connected, and 
w^^\i.) represents a loop operator corresponding to a loop created by the matrix 
A-^^ in the multi-matrix model. The operator ?^(o") describes an operator which 
changes the 'spin' on a loop locally from 1 to 2. Also w^^^ {£; [Hi^a)]^) describes a 
loop with [7i{a)y inserted. The symbol ~ means that as a function of £, the quantity 
has its support at ^ = 0. 

From eq. (|4.27|) for j = and n = 1, we have the relation 



+£i p di[ {w^'\i[)w^'\ii - ^^)w^'-\ll2)] 

+2giii2 {w^^\h + ^2))' ~ . (4.29) 



The planar part of the above relation agrees with eq. ( [4.18 ). Note that the loop 
amplitudes in eq. (|4.18| ) represent connected correlators. 



This agreement implies that Ti would correspond to B2- Taking into account the 
fact that Bn {n = Q mod m) do not exist and eq. ( |4.28| ), it is legitimate to consider 
that the amplitude (for j = 1, ■ ■ ■ , m) 

{w+{h)---w+{L)Bj) (4.30) 

corresponds to the connected part of the amplitude 

^ f da, («;«(^i) • • ■w^'\i,^^)w^'^ (£,; [nia,)y-') w^'\i,+,) ■ ■■w^'\in)) . (4.31) 



j=i 



5 Summary 



In this paper we have examined the correlators in the {m+l,m) unitary minimal 
models coupled to two-dimensional gravity from the point of view of the two-matrix 
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model. From the two-loop correlators and the wave function of the scaling operators, 
we derived the explicit form of the expansion of the loops in terms of the scaling 
operators. Using this expansion, we deduced the three-point functions from the 
three loop operators, and showed that simple fusion rules exist for all of the scaling 



operators. The three- loop correlator |13| can be understood to express these fusion 
rules in a compact form. 

At the (m + l,m) critical point in two- matrix models, the scaling operators Bj 
[j = mod m + 1) have no counterparts in the BRST cohomology of Liouville theory 



coupled to the corresponding conformal matter. In [|I^ ,these operators were argued 
to be boundary operators which couple to loops in the case of the one-matrix model. 
It was also shown explicitly that one of them, corresponding to a^+i in the case of 
the unitary matter, is a operator which measures the total length of the loops. 

We examined the role played by the rest of these operators. We showed, in some 
examples, that the operator Bn couples to the points to which n parts of the same 
kind of loops are stuck each other. In other words, the operator Bn connects n parts 
of the same kind of loops together. We think these operators play an important role 
concerning the touching of the macroscopic loops. The emergence of these operators 
in matrix models can easily be understood from the viewpoint of macroscopic loops 
and their expansion in terms of local operators. 
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Appendix A 

Let us prove that the two- loop correlator eq. ( |2.2|) can be written as eq. (pl7|). 
From eq. ( |2.2| ), we have 

"^-1 /sin TT— \ ^ 

w^h)wHi2)) = E (±)'"' —7^ {w,{i^)w,{e,)) (A.i) 



k=l 



7l/2 



and 



^ {Wk{£i)wk{i2)) = --^iii2 K,_^{M£,) K,^^{Me2) , (A.2) 



where we have introduced loop operators Wk{l) which represent loops with some 
distinct matter boundary condition. Making use of a formula 

KM KM-) = -J^ J K {-} exp (^-- —j (A.3) 
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and replacing t with tM^, we have 



1 /■°° dt 

2 Jo t ^"- \ t 



tM2 



exp 



??+£r 



2 2t ; 

Carrying out the integral with respect to M, and from eq. ( |A.2D , we have 

{Wk{il)Wk{i2)) 



1 f^ dt iii2 ,^ 

/ -fVi_jL 

Am Jo t t 






tW 



exp 



'i+^r 



2t 



Due to a formula, 



zK^_\^\{z)=j dE- 



i? sinh -kE 



10 cosh 7r_E — cos up 

the right hand side of eq. (|A.5|) turns into 

1 /"o^ dt 
Am Jo t 

Z7'„;„l 77" ^ f f ^ ' 



K,Az) 



/■°° , ^ E sinh vri? 

X / dE — - — 

JO cosh 7r£/ — cos vrp 



K. 



tM^ 



iE 



t 



exp 



2t 



Using a formula eq. (|A.3| ) again, eq. (|A.7|) turns out to be 



1 /■°° 
— / dE- 



E sinh vrii^ 



2m Jo cosh ttE — cos n— 



Putting eq. ( |A.8|) and eq. (|A.1|) together, we have 



-^(«-*(«)-si(±f-(^)^ 



A:=l 
/•oo 

X / dE 
Jo 



cosh IT E — COSTT — 



(A.4) 



(A.5) 



(A.6) 



(A.7) 



(A.8) 



(A.9) 



The integral in E can be carried out by deforming the contour. The residues for 
poles E = ±i{^ + 2n), n = 0, ±1, ±2, ■ ■ ■ , contribute to the integral and, after all, 
we obtain the following expansion for the two- loop correlators (for ii < £2) 



W+iii)w^ii2) 

1 m— 1 00 



£:E E (± 



>fc-i 



m 



k=l n=— 00 



k 

— + 2n 



l±+2n\iMil) Kj,^,jMi2) . (A.IO) 
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